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To Victor Petrovich Khavin 

Abstract. In this paper we give a simple proof of the fact that the average over all dyadic 
lattices of the dyadic //^-norm of a function gives an equivalent _ff^-norm. The proof we 
present works for both one-parameter and multi-parameter Hardy spaces. 

The results of such type are known; cf. [2] for the one-parameter case. Also, by duality, 
such results are equivalent to the "BMO from dyadic BMO" statements proved in [4| (one 
parameter case) and [S] (two-parameter case). 

While the paper generalizes these results to the multi-parameter setting, this is not its 
main goal. The purpose of the paper is to present an approach leading to a simple proof, 
which works in both one-parameter and multi-parameter cases. 

The main idea of treating square function as a Calderon-Zygmund operator is a common- 
place in harmonic analysis; the main observation, on which the paper is based, is that one 
can treat the random dyadic square function this way. After that, all is proved by using the 
standard and well-known results about Calderon-Zygmund operators in the Hilbert-space- 
valued setting. 

As an added bonus, we get a simple proof of the (equivalent by duality) inclusions BMO C 
BMOd, in the multi-parameter case. Note, that unlike the one-parameter case, 

the inclusions in the general situation are far from trivial. 
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Notation 

H^K, tensor product of Hilbert spaces % and /C; we consider its endowed with the 
"canonical" norm, which makes % ® IC & Hilbert space. If one of the spaces is 
a function space, for example if 7^ = L^, then 'H<^ K, can be interpreted as Lp' 
with values in /C. 

L^(X; %) L^-space of functions on X with values in %] X is usually with the Lebesgue 
measure. Symbols X and/or T-L can be omitted if they are clear from the context. 
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dx Lebesgue measure in if a; G 

V Dyadic Lattice in M^, see below. 

H^, real variable Hardy space on and its dyadic counterpart, associated to a 

dyadic lattice V. Wc will also consider Hilbert-space-valued spaces, and the 
notation like H^{X\T-L) will be used. 

^(Q)) \Q\ sidelength and volume of a cube Q C M^. 

expectation of a random variable ^, = j'^^(uj)dP(uj); sometimes, to distin- 
guish the random variable in a formula, the notation E^(w) or E^^{lo) will be 
used instead of E^. 



Cubes and dyadic lattices. Throughout the paper we will speak a lot about dyadic cubes 
and dyadic lattices, so let us first fix some terminology. A cube in M.^ is an object obtained 
from the standard cube [0, 1)-^ by dilations and shifts. 

For a cube Q we denote by i{Q) its size, i. e. the length of its side. Given a cube Q one 
can split it by dividing each side in halves into 2^ cubes Qk of size £{Q)/2: we will call such 
cubes Qk the children of Q. 

For a cube Q and A > we denote by XQ the cube Q dilated A times with respect to its 
center. 

Now, let us define the standard dyadic lattice Vq: for each k € Z let us consider the cube 
[0,2*^)-^ and all its shifts by elements of with coordinates of form j ■ 2^^, j G Z. The 
collection of all such cubes (union over all k) is called the standard dyadic lattice. 

A dyadic lattice V is just a shift of the standard dyadic lattice Vq. A collection of all cubes 
from a dyadic lattice X> of a fixed size 2*^ is called a layer of the lattice. 



Random dyadic lattice. Our random lattice will contain the dyadic cubes of standard 
size 2*^ (fc G Z), but will be "randomly shifted" with respect to the standard dyadic lattice 
Vq. The simplest idea would be to pick up a random variable x uniformly distributed over 
M.^ and to define the random lattice as x + Dq. This idea works for the torus T^, but 
unfortunately, there exists no such random variable x in M^, so we have to act in a little bit 
more sophisticated way. 

Let us construct a random lattice of dyadic intervals on the real line R, and then define a 
random lattice in R'^ as the product of the lattices of intervals. 

Let ill be some probability space and let x{lo) be a random variable uniformly distributed 
over the interval [0, 1). 

Let be random variables satisfying P{^j = +1} = P{^j = —1} = 1/2. Assume also 

that x{Ljj),^j(Lo), j G N are independent. Define the random lattice ^(u}) as follows: 

(1) We require that Io{u!) := [x{u) — l,x{u)] G T>{u!); this gives us all intervals in T){u!) 
of length 2'^, k<0. 

(2) To determine the rest of the intervals, it is enough to know dyadic intervals /^(a;) D 
Io{oj), of length 2^, k > 0. The intervals Ik{(^) are determined inductively: if 
Ik-i{co) G ^{oj) is already known (and thus all intervals of length 2^~^ in 'D{oj)), 
then 

• Ik{^) is the union of Ik-i{oj) and its right neighbor if ^fc(i^) = +1, and 

• Ik{^) is the union of Ik-i{oo) and its left neighbor if ^fc(^) = — 1- 

To get a random dyadic lattice in we just take N independent random dyadic lattices 
Pi, • • • , Vn in R and consider all cubes Q = Ii x I2 x . . . x In, G V^. 
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0. Introduction and main results 

0.1. One parameter case. Let = i?^(M^) be the usual real variable Hardy space on 
M^, and let be its dyadic counterpart, defined as follows. 

Consider a dyadic lattice T> in M^. Let = E^ be the averaging operator over cubes 
Q € P of size 2^^, Ekf{x) = \Q\"^ Jq f, where Q is the cube in V of size 2^ containing x. If 
QeVisa cube of size 2^= define E^ by E^f = {\Q\-^ f)!^ = lQ{Ekf). 

Define the martingale differences = A|' := E^_^ — E'^, and again for a cube Q ^ T> of 
size 2^= define by A^/ = lQ{Ekf). 

Define the dyadic square function S = by 



(0.1) (5^/)(x)= fj]|AF/(x)|') =1 



1/2 



Vfcez / \QeV,QBx 

One can also consider a slightly different square function Sf = S^f 

1/2 

7V I aI>^|2 



1/2 



(0.2) (5^/)(x)= j;(i?|'|AF/|^)(.) = [EQKf\ 

VfceZ / \Q(^V,Q9x 

in other words, each term in the sum on the right equals the average of |Aq/P if x G Q and 
otherwise. _ 

It is not hard to show and will be explained later that C^^||5/||i < US'/Ill < C||5/||i, 
where the constant C depends only on the dimension N. 

Definition. A function / € Lj^^ is in the dyadic Hardy space (with respect to the dyadic 
lattice V) if il/ll^i^ := < oo 

One can use the square function S-rif here and get an equivalent norm. 

And now one of the main results of the paper. Let I'(w), a; G Q be the random dyadic 
lattice, as described above, and let us recall that E = E^^ denotes the expectation (average 
with respect to co) 

Theorem 0.1. A function f € ^[^^^(M") belongs to if and only if 

[mSv(..)fix)\^)Y^^dx<oo. 

Moreover, the latter quantity gives an equivalent norm on . 

The same result is true with the square function Sx>[uj) in place of Sx>{uj) ■ 

0.2. Multi-parameter case. The above results can be generalized to the case of multi- 
parameter Hardy spaces. Let Xj^ = M^*^, A; = 1, 2, . . . , n and let H^{Xi x X2 x • • • x Xn) be 
the n-parameter Hardy space, see Section [2.11 below for the definition. 

Define the dyadic Hardy space as follows. Let Pfc be a dyadic lattice on X^, k = 1,2, . . . ,n 
and let V = Vi X V2 X . . . X Vn he the product dyadic lattice on X = Xi x X2 x . . . x X„; 
the elements on V are the "rectangles" (parallelepipeds) R = Qi x Q2 x . . . x Qn, Qk £ T^k- 

For a multiindex k = {ki,k2, ■ ■ ■ kn) define on X = Xi x X2 x ••• x X„ the average 
Ek '■= E\^E1^ . . . E^^ and the martingale differences A^ := A^^ A|^ . . . A^^, where E^^, and 

A^. are the "one variable" averages E^^ martingale difference A^^. as defined above by in 
Section 10. H taken in the variable Xj € Xj . 
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For a "rectangle" R = Qi x Q2 x . . . x Qn, Qk ^ T^k define := Aq^Aq^ . . . Ag^, where 
again Aq, is the operator Aq. defined by (j0.2p taken in the variable Xj. 
Define the multi-parameter square function S = Sx> by 



\ 1/2 

(0.3) (5„/)(x) = I > : |Ak/(x)p = I > : \ARf{x)\' 




1/2 



|2 



1/2 

One can also define the square function S = Sx> 

(0.4) Sf{x):=(Y.{^k\Akf\'){x)\ =[ ^ {ER\ARf\ 

The definitions look very similar to (|0.1|) . (|U.2|) . only here the sums in the right hand side is 
taken over all "rectangles" R, while in (|U.1|) . (|U.2|) they are taken over all cubes. 

We use the same notation for the one-parameter and multi-parameter square function, but 
since we will treat these cases in different sections, we hope to avoid the confusion. 

Definition. Let V = 'D1XV2X . . .x P„ be a product dyadic lattice on X = Xi x X2 x . . . x X„. 
We say that a function / G LI^^{X) belongs to the dyadic Hardy space if := 
\\Sj,f\\i < 00 

Let now 'Di{uj),'D2i^), ■ ■ ■ , T^n{'^) be the independent random dyadic lattices on the spaces 
Xi,X2,...,X„ respectively, and let T){uj) = T)i{oj) x T>2{oj) x ... x T>n{oo) be the multi- 
parameter random dyadic lattice. 

Theorem 0.2. A function f G ^[^^(M") belongs to H^{Xi x X2 x ■ ■ ■ x X„) if and only if 

f [n\Svi.)fix)\')f'dx<oo; 

here Sx)(ij) is the multi-parameter random dyadic lattice defined above. Moreover, the latter 
quantity gives an equivalent norm on . 

The same result is true with the square function Sxi[uj) in place of Sxi[uj) ■ 

0.3. Some remarks. 

Remark 0.3. It is a simple exercise to show that in both one-parameter and multi-parameter 
cases 



n\Sv{u)f{x)\^ = / \Svf{x + u 



)\^du 



where P is a fixed dyadic lattice and the limits is taken over cubes (or parallelipipeds) Q C 
{Q <Z X ui the multi-parameter case) centered at whose sidelengths tend to 00. 

Remark 0.4. Note that for every dyadic lattice V one has < ells'!)/!!!, thus by Tonelli 

theorem and Holder inequality we have for p <2 



< C I [E(|5p(^)/(x)r)]'/^dx < C / [E{\Sr>^^)f{x)\^)f'dx. 



In the one-parameter case the estimate ||/||//i < CHS*©/!!! is trivial and well known. Indeed, 
H^-BMO duality {{H^)* = BMO, (H^)* = BMOd, see the definitions below in Section 
and the trivial inclusion BMO C BMOx) imply the C with the desired estimates of 
the norms. 
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Since BMO 7^ BMO^, one can conclude that th inclusion C is proper. 

In the multi-parameter case the estimate < C||5x)/||iis also known to specialists, 

but it is less trivial. In fact, the only place that the author is aware of, where this is proved 
is the Ph.D. thesis of J. Pipher; this proof is far from trivial and the calculations are quite 
tedious. Below in Section 12.31 we present a different, quite simple prove of this fact. This 
proof is based on one-parameter (Hilbert-space-valued) //^-BMO theory. 

Remark 0.5. Applying Holder inequality to Theorem 10.11 (to Theorem 10.21 in the multi- 
parameter case) we get that for p <2 

(0.5) / [E(|5p(,)/(x)r)]'/^dx< / [E(|5^(,)/(x)|2)]^/'dx<C7||/||^, 

If p = 1 Tonelli Theorem and (|0.5p implies that 

(0.6) E||5^(^^/||i<C||/||^i. 

This was proved in one-parameter case in [2j. Note, that by duality (|0.6p is equivalent to 
"BMO from dyadic BMO" statements: if /(^ w G O is a measurable family of functions, 
GBMOp(^), 11/^11 

bmOt,, — then for / defined by f{x) — Etj/i^(x) one has / G BMO, 

II/IIbmo < c. 

This "BMO from dyadic BMO" result was proved directly in [3] (one-parameter case) and 
in the recent paper ^ (the two-parameter case). 

1. Proof of Theorem 10.11 

1.1. "Vectorization" of the square function. There is a standard way of making the 
nonlinear operator Sx> into a linear one by treating Sx>f as a vector-valued function. 
Namely, let us consider the space L^(£^) consisting of functions on Z x such that 

/ \ 1/2 

\\f\\LHi^)--= [ {T.\fik,x)n dx<c^. 

■'^'^ VfcGZ / 

Define the vector- valued square function Sj?/ by 

Svf{k,x) := Aff{x), k£Z,x£R^ 

{Svf is a function on Z x M^). We will treat this function as a function of the argument 
X G with values in f = 

Clearly, | Sx>/( • , a;) | ^2 = Sx>f{x), so / G if and only if S/ G L^{£'^). Moreover, 
I|Si?/||li(£2) = \\Svf\\i = ll/IIz/i,. 

Let now I?(w) be the random dyadic lattice, and let 0,P be the corresponding probability 
space. Consider the space C = L^{i'^ (g) L^(0,P)), 

\ 1/2 




Y,\fik,uj,x)\^dFiuj)] dx. 



^ feez 



It is an space with values in the Hilbert space l"^ ® L^(il,P); here again / is treated as 
function of the argument x G with values in (8> L^(0,P). 

Define the vector- valued square function S with values in (8) L^(ri,P) by 

(1.1) Sf{k,uj,x) = ST,(u;)fik,x), X eR^, k eZ,u; en. 

here and below we will use notation S( . , ■ ,x) =: S/(x) G (X> L^(il, P). 
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Clearly, \\Sf{x)\\(,2^L2(^^^p) = [E{\Sx>(oj)f{x)\'^)Y^^,so ||S/||^ = Jj^jv [E([5'x)(^)/(x)p)] dx 
(the norm in jC was constructed so this would hold). 

Similar "vectorization" can be performed to the square function S. Namely, fix some 
ordering of the "children" of a dyadic cube (the same one for all cubes), and for each dyadic 
cube Q define operator from the functions constant on the "children" of Q to C by 

UqIq, = \Qk\^/^ek, k = l,2,...,2^'', 

where Qk are "children" of Q and {efc}|^^ is the standard basis in M^^. Let now £ := 

L'^{f ® L'^{n,¥) (g) C^"^) with the norm 

1/2 



■= [ if ^\mu;,x)\\^dF{u)] dx. 
Define the square function S by 



where Q € '^{oj) is the cube of size 2^ containing x. From the construction it is clear that 

||S/||^ = V[E(|^^(,)/(x)p)]'^'dx. 

Now the proof of Theorem 10.11 can be outlined in few sentences. First, it is nor hard to 
show that S (or S) is a Calderon-Zygmund operator, whose kernel takes values in the Hilbert 
space £^ (8) L^(i7, P). It is a well known fact that such C alder 6n-Zygmund operators map 
to L^, so 

f [n\Sv(.)f{x)?)]"^dx=\\Sf\\^<C\\f\\H^. 

Jrn ^ 

As we discussed above in Section \07l] the opposite inequality < C||S/j|^ is trivial. 

Of course, the classical theory of Calderon-Zygmund operators deals with the scalar- valued 
kernels. But, as it is well known to the specialists, all the facts that we need, are valid in the 
case of Hilbert space valued kernels too. 

However, the blind trust is not expected from the reader: all relevant facts will be presented 
below. 



1.2. S as a Calderon-Zygmund operator. Let us recall that the classical Calderon- 
Zygmund kernel on is the function A'( • , • ) defined on X \ {{x,x) : x G R^} 
satisfying 

(1) \K{x,y)\<C\x-y\-^; 

(2) There exists 6 > such that 

\K{x, y) - K{xo, y)\, \ K{y, x) - K{y, xo) | < C^r ~ 



xo\^+^ 



whenever \y — xq\ > 2\x — xq\ 

One can also consider operator-valued kernels, K{x,y) € B{X,Y) for arbitrary Banach 
spaces X and Y. In this case | • | in the left hand side should be replaced by the norm in 
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1.2.1. Vector-valued Calderon-Zygmund operators. A Calderon-Zygmund operator with (op- 
erator-valued) Calderon-Zygmund kernel is a bounded operator T : L^{R^;X) 
L^{R^,y) such that for all compactly supported / G L'^{R^;X) and g G L'^(R^;y*) with 
separated supports 

{Tf,g)= [[ {K{x,y)f{y),g{y))dydx. 

Such operators with operator- valued kernels were considered, for example, in [9], and it was 
proved there (see Theorem 1.2 in Ch. 3) that if X and y are Hilbert spaces then the operator 
acts from to L\ WT/W^i^kn^p^) < C\\f\\Hi(^Nyy 

In fact, in [9] a much more general situation was considered. The kernel was only assumed 
to satisfy a weaker version of condition (2), and there was no condition (1). Moreover, it was 
assumed that T was bounded in some LP, 1 < p < oo, not necessarily p = 2. 

Note, that condition (1) in some form is required for the proof of Tl and Th theorems, but 
if one assumes that T is bounded in some L^, 1 < p < oo the condition (2) alone is sufficient 
for the action from to L} . 

Also, in ^ the theorem was proved for the case when X and y are arbitrary Banach 
spaces, if one defines via atomic decomposition. 

It is well known that for the case of Hilbert-space- valued functions all the definitions of 
(via atomic decompositions, via maximal function, via different square functions, via Riesz 
transforms) are equivalenlQ, so one can use the result from [9] without worrying about what 
definition of is used. 

In this paper we are considering the case when X = C and y = H = i'^ L^(J7, P), so we 
can say that K takes values in the Hilbert space H. 

Operators with such kernels also act naturally from L^(R^; /C) — )• L^(M^; 7^ (8) /C), where 
/C is a Hilbert space, and we will need this interpretation later in Section [2j Indeed, with 
each vector h € Ti we can associate an operator K.^f^f(^h^K,®'H (and the norm 
of this operator is \\h\\). So, if K is ah ?^-valued Calderon-Zygmund kernel, then a Calder- 
on-Zygmund operator on L2(M^; /C) is a bounded operator T : L'^{R^; K) L'^{R^; U^IC) 
such that 

{Tf, g) L\R« ■,'H^K.) = II {K{x, y) (S) f{y),g{y))y,(^K: dydx. 

for all compactly supported / and g with separated supports. 

1.2.2. Why S is a Calderon-Zygmund operator? To find the kernel X of S we need to compute 
S6y, where 6y is the unit mass at y G R^: 

K{x,y) = \A^^''^6yix)] ee®L\n,¥) 

To get that expression rigorously, one needs to approximate 6y by appropriate bump func- 
tions; notice that /S.^'^dy is well defined, i.e. does not depend on the choice of approximating 
sequence. 

Notice that for a dyadic cube Q, IS.Q5y{x) / only if x G Q, y G Q. Therefore Aq5j,(x) = 
if i{Q) <\x- y|oo, so Al^'"^5y{x) = if 2^= < |x - y\^. 

Noticing that ||A^"<5y||oo < C2~^^ and summing the geometric series we get the property 
(1) of Calderon-Zygmund kernels. 



Unfortunately, the author cannot point to a paper where all such equivalences are proved; but following 
the proofs for the scalar-valued case, one can see that everything works for the case of Hilbert-space valued 
case as well. 
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To show property (2) notice that A^^'^^ 6y{x) = A^^'^^ 6y{xo) if all 3 points x,xo,y are in 
the same cube Q G ^{uj), £{Q) = 2^ and the points x, xq are in the same "child" of Q. 
The probability that it fails for a given k can be estimated above by C\x — X(}\2~^. Using 



the estimate ||A^"5y||oo < C2 we conclude that 



/ lA^^'^hyix) - A^'^''''6y{xo)\^dF{uj) < C2-'^^^\x - xo|2~^ = C\x - xq\2- 
Jn 



-2kN~k 



Let kf) be the maximal A; € Z such that 2^ < minjlx — y|oo5|2;o — y|oo}- Then A'^'^6y{x) 



Af'^6y{xo) = for /c < ko, so 



\K{x,y)-K{xo,y)\' < C 2-''=^-'lx - xqI < C|x - xo|2-2^'«^-^o 



-2N-1 



k>ko 

< C\x - xo\\xo - y\ 

Interchanging x and y and repeating the above reasoning we also get that 

I K{y, x) - K{y, xq) | ' < C\x - xo\\xo - y\~^''-\ 

This means condition (2) holds with 6 = 1/2. 
The proof for S is absolutely the same. 

1.3. A remark about conditions SI = 0, S*l = 0. Material in this section is not needed 
for the proof of the main results. However, it might be of interest for specialists; one can 
use it to present a different proof of the main results, without employing the cited above in 
Section fl . 2 . 1 1 r esult from [9j result about Calderon-Zygmund operators with operator- valued 
kernels. 

Note that the operator S, introduced above satisfies the conditions SI = and S*l = 
(more precisely, the second condition should read as S*le = for all e G "H), which are well 
known to everybody familiar with T(l)-theorem. 

If one formally plugs 1 into S or S*, the result will be 0. Of course, it is only a formal 
reasoning, for 1 is not in the domain of S, but any reasonable interpretation of SI gives the 
same result. For example it is not hard to show that 

(1.2) SIq^O, Sn^e^O asi{Q)^oo 

uniformly on compact subsets, where cubes Q are centered at 0. 

It is also easy to see that SI = 0, but unfortunately S*l ^ 0. However, it is easy to modify 
S to make S*l = 0. 

Namely, let if he a function on the cube [0, 1)^ taking values ±1 and such that Jg ipdx = 
and let "^Qix) = ^{{x — Xq)/^{Q)), where x^ is the base of Q, i.e. the point in Q with 
smallest coordinates. 

Define the square function S by 

Sf{k,uj,x) :=<^g(x)C/QAQ/GC2 , 

where Q G 2^(w) is the cube of size 2^ containing x. 

The function ip^ in the definition of S is introduced to insure that S*l = 0. Now it is easy 

to show that SI = 0, S*l = (in the sense of ([LSI))- 

Calderon-Zygmund operators satisfying Tl = and T*l = map H'^ H^. To show 
that one, for example can consider matrix of such an operator in the wavelet basis. It was 
shown in j6] that under rather mild assumption about wavelet basis, the coefficient space of 



AND DYADIC 



9 



in this basis is the TViebel-Lizorkin space f^' ; see Appendix (Section [3|) for the definition. 
In [6] the scalar-valued case was treated, but one can easily see that everything works for the 
Hilbert-space valued case, and one just get the vector-valued space f^' . Moreover, in |5j the 

spaces with values in UMD Banach spaces were characterized in terms of coefficient in 

'02 

the wavelet basis; in the Hilbert-space- valued case the result gives exactly f^' with values in 
the Hibert space. 

Using the standard estimates with Calderon-Zygmund kernels one can see that if a Cald- 
eron-Zygmund operator T (even with the operator- valued kernel) satisfies Tl = 0, T*l = 0, 
then its matrix in the wavelet basis is what is called in [3j almost diagonal for fp''', \ <p,q < 
oo; cf. Section [3] below for the definition. 

And it was shown in [3j that almost diagonal operators are bounded on all fp'^, 1 < p,q < 
oo. Since the almost diagonality is a condition on the magnitude of the entries, the result 
holds for vector-valued Triebel-Lizorkin spaces. Of course, instead of considering a wavelet 
basis, one can consider a frame decomposition, given by what is called in [3J ip -transform; all 
the estimates will be the same. 



2. Proof for the multi-parameter case 

Proof of the main result for multi-parameter case (Theorem 10. 2p follows the lines of the 
proof for the one parameter. 

Without loss of generality we can assume that the probability space 17 is represented as a 
product, (r2,P) = (r^i X X • • • X ^nilPi X P2 X ... X P„) and that the random dyadic grid 
T>k{oj) on Xk depends only on ojk- 

For a dyadic lattice P = x D2 x . . . x D„ on Xi x X2 x . . . x X„ define the vector-valued 
square function S-p taking values in l'^{'L^) 

So/(k, x) := A'^fix), k G Z", X G 

Consider the space C = L1(£2(Z") L^{n,F)), 

\ 1/2 

Y,\fiKi^,x)\^dF{u;)\ dx. 
^ keZ" / 

It is an space with values in the Hilbert space £^(Z") (g) L^(r2,P). Note, that this Hilbert 
space can be decomposed as ^^(z") P) = (i^ L2(17i, P^)) (^2 ^ l2(J72, P2)) • • • O 
L^nn,^n)y, here f = fiZ). 

Define the vector-valued square function, taking values in the space ^^(Z") L'^{Q,¥) by 

Sf{k,Lo,x) = S2,(^)/(k,x), xeR^, k G Z",a; G 17. 

Note, that S can be decomposed as a tensor product S = Si (gi S2 (X" • • • (81 S„, where Sfc is the 
one parameter square function defined by (jl.ip in variables Xk G Xk = M^*, oJk € ^k- 

Clearly, as in the one parameter case, we have for the multi-parameter square functions 
IIS/Ill = [E^(|5^(^)/(x)|2)]i/2 dx. ^ 

Similarly, for the square function S one can construct its vector version S with values 

in n = Til ^ n2 ... ^ Un, Tik = i'^ ^ L'^{^lk,'^k) ^ C2'^\ Again, it is easy to see that 

1/2 




- fx 



IE.(|52,H/(x)h 



dx. 



As it was already discussed in Section 11.21 operators Sk are (Hilbert-space- valued) one- 
parameter Calderon-Zygmund operators, so S is the tensor product of such operators. And 
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it is probably immediately clear to experts, that such operators map H^{Xi (g)X2 (8 • • • 
to L\ 

One way to see that is to notice that S is a trivial case of multi-parameter Calderon- 
Zygmund operators, and according to Theorem 2.2 in |7j such operators map H^{Xi (g) (8) 
. . . (8) Xn) to L^. Of course, one needs to use a Hilbert space valued version of the theorem, 
but it is clear to the specialists, that the proof from [7] works in this case. It is also clear that 
while Theorem 2.2 in [7] is stated for M x M x . . . x M, the proof works for (g) -^^2 (8 . . . Xn, 
Xk = R^K 

For the reader who is not well familiar with multi-parameter spaces we present below 
an alternative proof, which exploits the tensor product structure of S = Si S2 . . . (g S„; 
it uses only theory of one parameter ff^-spaces. Of course, we will need the theory of H^- 
spaces with values in a Hilbert space, but we need the vector valued theory in the above 
multi-parameter reasoning as well. 

There is one more reason for the presenting the one-parameter proof below: the reasoning 

1 /2 

above gives the estimate [Ki^{\Sx>(u})fix)\'^)] dx < C\\f\\Hi(^x)- The opposite estimate 
follows from the inequality < \Svf{x)\dx, which is trivial in one-parameter case. In 

multi-parameter case, the same estimate, while true and known to specialists, requires some 
some work to prove it. The one-parameter approach presented below gives a reasonably 
simple proof of this estimate. 

2.1. Multi-parameter i7^-spaces. Recall, cf [TOl III. 4. 4] that for for one-parameter Hardy 
space H^{M.^) the norm ||5'^/||i, where is the Lusin square function (aka Lusin Area 
Integral) 

5V(x)= / |V/(y,t)|V-"dydt 

(here Tx '■= {{y,t) : y € M^,t > 0, |y — j;| < t}, and f{y,t) is the harmonic extension of / 
from to R^+^ := x M+ = {{x,t) : x G R^,t £ M+}) gives an equivalent norm. 
One can consider the vectorization of as follows. Let F = Fq and define 

SV(^, y, t) = t(i-")/2V/(x + y, t), X G M^, (y, t) G F. 

By construction S^/(x, •, •) G -^^(F) and ||S'^/(x, •, ■)\\L2(r)®c^ = S^fi^)^ therefore 
l|SV||Li(L2(r)) = II^Vlli. 

One can define multi-parameter square functions and by 

S^f {xi,X2, ...,Xn): = 

/ IV1V2... V„/(yi,ti,y2,t2,---,yn,in)px 

nl/2 

X tl'^H^^"^^'' ■ ■ ■ ti''^"dyidtidy2dt2 . . . dyndtn ; 

here /(yi, ti, 2/2) ^2) • • • 1 ym ^n) is the harmonic in each variable {yk,tk), yu £ M^'', G M4. 
extension of / from M^i x M^^ x . . . x M^" to M^^^^ x M^^+i x . . . x M^"+i and is the 
gradient in the variable {yk,tk)- 

Following [1] we say that / G H^{X) = H^{Xi ® X2 ® . . . ^ X„) if ^/ G L^{X) and 
||S^/||i defines one of the possible equivalent norms in H^{X). 



AND DYADIC 
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We also define the vector- valued linear square function as Sf{x) = ® S2 <8> • • • ® 
S^/(x) £ n = Til 0712(8) ■■■ (S)n„, where it one-parameter square function defined above 
taken in the variable x^, Tik = L'^^Xk) and is the cone in M^*"^^ with the vertex at 0. 

Again, by the construction |[S^/(x)||-^ = \S^f{x)\, so ||S'"/||/,i(x.-^) gives the norm in 
H\X). 

We will use to define multi-parameter with values in a Hilbert space /C; in this case 
€ ^ X /C. We introduced such spaces only for notational purposes, so while most of 
the theory of multi-parameter //^-spaces can be transfered to the Hilbert-space-valued case, 
we do not need this. 

1 /2 

2.2. Proof of estimate/^ [E^(|52,(^)/(:e)|2)] dx < m{x)- Consider a multi-param- 
eter square function S2 C?) • • • (8> S^, where each S^, 2<A;<nis either one-parameter or 
one-parameter "random" square function S, defined in (|l.ip . taken in the variable x^. 

Assume that the choice of is fixed. For a scalar-valued / the function 82 ® ■ ■ ■ ® Snf 
takes values in := <S) ■ ■ ■ fS) fin, where each fik is either L'^(Tk) or i'^ (g) ^^(1^^, P^), 
depending on what square function is. 

Let S^, S/c be the Lusin and "random" square functions, taken in the variable x^, and let 
Tik ■= L^(Fjt) (g) and Ti'j^ := i'^ (gi L^(r2fc,Pfc) be the corresponding target spaces. 

Lemma 2.1. 

|Si«)S2«)...®S„/(a;)|| dx<C i|S^^ «) §2 ® . . . (g S„/(x)|| dx 

Since the tensor products of square functions we consider does not depend on the order 
(the square functions, taken in different variables obviously commute), the above lemma tells 
us that one can replace a factor by Sk in Si (g) S2 fg) • • • fg) S„/ and increase the norm by 
at most the factor C. 

Starting with (g . . . (g S^/ and applying Lemma [2?T] successively to each factor, we 
get 

/ ||Si (g S2 . . . (g Snf{x)\\^,dx <C I S^f{x)\\^dx, 
which is exactly the desired estimate (here H' = T-L'i (g) 7^2 "SD • • • (X" V-'n and % = Hi ® 7^2 ® 

Proof of Lemma \2.1\ Let us introduce notation x = (2:1,3;-'^) € X = Xi x X"^, where x"^ = 
(X2, ^3, . . . , x„) G := X2 X . . . X X„. 
Consider the vector-valued function 

Sj^ ® S2 (g . . . ® s„/ = sl^ O (S2 ® . . . ® S„)/. 

If 



/ \\S\(S){S2(S)...®Sn)f{x)\\ dx< 



00, 



we conclude that for almost all x^ 

(2.1) "^2® . . . (8)Snf{ - ,x^) £ H\Xi-H^), 
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and 

(2.2) / ||SL0(S2 0...®S„)/(x)|| 

Jx Hi<»H' 

=/^jrs,s...»s„/(.,.')ii^.^^_^^.^</.^ 

Note, that we have in (|2.ip the usual one-parameter //^-space (although vector- valued) . 

As we discussed above in Section [T72[ Si is a vector- valued Calderon-Zygmund operator, 
so it maps one-parameter to (even in the Hilbert-space- valued case), so for almost all 

we have 

(2.3) [ ||Si0(S2®...^S„)/(xi,xi)|L , -^Jxi 

JXi H^^H 

<C||S2^...0S„/(.,xi)||^^^^^^^^^ 

= C f \\S\®{S2®...®Sn)f{xi,X^)\\ --^dxi 
Jxi Hi<»W 

Integrating over and taking into account (j2.2p . we get the conclusion of the lemma. □ 

2.3. Estimate < C* \SDf{x)\dx. Proof of this estimate follows the lines of Sec- 

tion 12.21 almost word by word. It is based on the following analogue of Lemma 12. H which 
allows us replace one-parameter Lusin square functions by the dyadic ones. 

Consider again a multi-parameter square function S2 ^ . . . (8i S^, where now each S^, 
2 < A; < n is either one-parameter or one-parameter dyadic square function Sx>^, defined 
in (jl.ip . taken in the variable x^- We assume here that in each Xf. dyadic lattices Vf. are 
fixed. 

We will use the same notation as in Section 12.21 with the only exception that now = i'^ 
Lemma 2.2. 

/ ||S^^«)S2®...®S„/(x)|| ~,dx<C I ||Sp, (»S2®...(8)S„/(x)|| ~,dx. 

Jx Hi®H Jx H-^®H 

Applying Lemma 12.11 successively to each variable, as we did in Section 12.21 we get 

j S\f{x)\\^dx <C j ||S©, Si5„/(x)||^,dx, 

which is exactly what we need. 

Proof of Lemma \2.SX Similarly to (j2.2p we get 

(2.4) / \\Sv,®{^2®...®Sn)f{x)\\ ~JX 
Jx V.i®H' 

= [ ||S2^...$5S„/(.,xi)|| ,^ ~,dx\ 

We will now use the fact that for Hilbert-space- valued functions Hv^Hi^i < C||(/3||j|^i^. Again, 
as in the scalar-valued case, it follows from H^-BMO duality {{H^)* = BMO, (H^)* = 
BMOx)) and the trivial inclusion BMO C BMOx), which imply the inclusion H.^ C with 
the desired estimates of the norms. 

Using this inequality we get that for almost all x^ 

1IS2 ® . . . ^ s„/( • , ^')\\^,^x,m ^ c'iis2 . . . ® s„/( • , 11^, • 



AND DYADIC 

Integrating over and using (j2.2p . (|2.4p we get the conclusion of the lemma. 
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3. Appendix: some facts about and BMO spaces. 



3.1. Hilbert-space- valued BMO spaces. Let us recall that a function on X 
values in a Hilbert space H belongs to the space BMO = BMO{X,'H) if 



pN 



with 



(3.1) 



BMO := sup / ||/(x) - /„|L<ia; < oo; 
Q Jo " 



here := \Q\ ^ fg f{x)dx and the supremum is taken over all cubes Q C 

If we fix a dyadic lattice "D and take the supremum in ()3.ip only over dyadic cubes Q (zV, 

we get the dyadic space BMOx) associated with this lattice. 

It is weh known that (F^(M^;H))* = BMO{R^;H) and {H^{R^ ■,n))* = BMOi?(M^; ?^); 

any standard proof of H^-BMO duality would work for the Hilbert-space- valued functions. 



3.2. Triebel Lizorkin spaces ip''^ and equivalence of HS/Hi and US'/Hi. In this section 
we fix a dyadic lattice P, for example take for V the standard dyadic lattice. 



3.2.1. Spaces fp . Following [3] define the space fp , (a G M, 1 < p,p < oo, consisting of 
sequences s = {sqIq^p such that 



1/9 



< oo 



LP 



For p = oo the norm is defined using BMO-like norm 



:= sup [ (iQr"^"i.^Qi • iqi-'^'^qA 

We are interested in the case when the smoothness parameter q = 0; to simplify the 

0,9 



notation in this case we will use fp := fp 



We will need the following facts about duality for spaces fp '"^i 
here 1/p + 1/p' = 1, 1/g + l/q' = 1. 



3.2.2. Almost diagonal operators. Following [3] we say that an operator A with matrix 
^^Q P^Q PeV almost diagonal (for spaces fp 
such that 



fp''') if there exists e > and C < oo, 



K,p\ < C I 1 + niax{£(P),£(Q)} 



-N-e 



X mm 



Km) 



{N+e)/2 



(N+e)/2 



(the definition is a bit more complicated for fp''' with s ^ 0) 

It was shown it [3j that an almost diagonal operator is bounded in all fp spaces, 1 <,q,p < 
oo. 
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3.2.3. Equivalence of \\Sf\\i and \\Sf\\i. From the above result one can easily obtain the 
equivalence of \\Sf\\i and US'/Hi- First, since maxQ |Aq/P < 2^ EgdAgfl"^) we have point- 
wise estimate Sf{x) < 2^/^Sf{x) and so ||S/||i < 2^/^\\Sf\\i. 

To get the estimate HtS/Hi < CyS"/ ||i let us express the conditions Sf £ L^, Sf E in 
terms of Triebel-Lizorkin space f(. Namely, with each function / G L^^^ let us associate 2 
sequences a = {a^l^^^ and b = {b^j^^^ 

1/2 



„\ -1 l/Z 



where R is the "parent" of Q and x is an arbitrary point in Q (the result does not depend 
on x). Then clearly 

ll^/lli = ll«llf^, l|5/lli = ll%- 

Note that 

\ Q is child oi R J Q is child of R 

where |6| := {{bQ^Q^^ and for s = {sq}q^^ 

iTs)^:=2--/^ E ^Q- 

Q is child of R 

The operator T is almost diagonal (it has only finitely many "diagonals"), so 
\\Sfh = Mi^^ < \\T\b\h^ < c\\ \b\ llf, = C\\b\\i, = C\\Sf\\,. 

□ 
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